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Abstract. This paper investigates the nature of the MHD turbulence excited by the streaming of accelerated 
cosmic rays in a shock wave precursor. The two recognised regimes (non-resonant and resonant) of the streaming 
instability are taken into account. We show that the non-resonant instability is very efficient and saturates through 
a balance between its growth and non-linear transfer. The cosmic-ray resonant instability then takes over and is 
quenched by advection through the shock. The level of turbulence is determined by the non-resonant regime if the 
shock velocity Vsh is larger than a few times ^cr c, where ^cr is the ratio of the cosmic-ray pressure to the shock 
kinetic energy. The instability determines the dependence of the spectrum with respect to fcy (wavenumbers along 
the shock normal). The transverse cascade of Alfven waves simultaneously determines the dependence in k±. We 
also study the redistribution of turbulent energy between forward and backward waves, which occurs through the 
interaction of two Alfven and one slow magneto-sonic wave. Eventually the spectra at the longest wavelengths 
are found almost proportional to k'^^. Downstream, anisotropy is further enhanced through the compression at 
shock crossing. 

Key words. Physical processes: Acceleration of particles - Shock waves - Turbulence - Interstellar medium: 
Supernova remnants 



1. Introduction 

Fermi acceleration of cosmic rays in astrophysical shock 
fronts depends in a crucial way on their transport prop- 
erties in the turbulent magnetic field on both sides of the 
shock. Often the turbulent field spectrum and intensity 
are arbitrarily prescribed, assuming that it has been built 
by the ambient medium i ndepend e ntly o f the shock accel- 
eration process. However iLerchd ()l967j) . IWentzell ()l969|) 
have argued for a long time that the development of an 
anisotropy of the cosmic ray distribution function triggers 
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an instability upstream of the shock. iMcKenzie fc Vol^ 
(|l982) have investigated the consequences of this phe- 
nomenon in the energy budget of the shock, in particular, 
with respect to the efficiencies of conversion of the kinetic 
energy into therm al, turbulent m agnetic and cosmic ray 
energies. Recently, iBell fc Lucekl l|200l() have shown that 
the amplification of the turbulent magnetic energy could 
be quite significant, producing a magnetic field intensity 
suitable to push the high energy cut-off of the proton dis- 
tribution up to the "knee" of the cosmic ray spectrum 
{Er^3x IQi^eV). 

This theory has then been developed further with accurate 
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inv estigations of supernovae rem nants (SNR). In particu- 
lar IPtuskin fc Zirakshvilil l)2003l) have analysed the gen- 
eration of turbulence and emphasized the importance of 
advection for the saturation of the spectrum. They have 
also carried out a preliminary examination of the role 
of the Kolmogorov cascade in th e ene r gy tra nsfer among 
the excited waves. More recentlv iBelll ()2004l) described a 
non-resonant regime of the streaming instability and has 
shown that its growth rate should be dominant in the high 
wavenumber range (to be discussed in more detail below). 
The fast growth of non-resonant modes could provide the 
necessary magnetic field intensity during the early stages 
of the SNR evolution to accelerate cosmic rays even up 
to the cosmic-ray spectrum "ankle" at £' ~ 3 x 10^^ eV. 
If verified this possibility would bring a strong support to 
th e standard galactic cosmic-ray model [see the discussion 
in iDrurv et al.l l)200lj) and references therein] . 

In this paper, we analyze the excitation of Alfven 
waves as a function of the location in the upstream flow 
and of the wavenumber taking into account the two insta- 
bility regimes (Section 2). In Sections 3 and 4, we calculate 
the saturation mechanism of the instability considering 
the advection effect as a function of the wavenumber and 
the location in the upstream flow. We calculate the contri- 
butions of two non-linear effects: the transverse non-linear 
transfer among turbulent Alfven waves, and the non-linear 
backscattering of Alfven waves off slow magneto-sonic 
waves. These two processes are shown to be relevant and 
essential to the determination of the anisotropic turbu- 
lence spectra. We finally derive these spectra which are es- 
sential to calculate the cosmic-rays transport coefficients. 
The detailed calculation of these trans port coefficients is 
carrie d out in the companion paper ijMarcowith et alV 
l2nn(i^ , hereafter paper II. In section 5, we examine the 
consequences in the downstream flow of the upstream ex- 
citation of the turbulence. In particular we propose a spec- 
trum for the turbulent magnetic field and estimates of the 
relaxation length and of the parameters describing the dy- 
namo action downstream. The technical derivations are 
presented in two appendices. 



2. Upstream excitation of MHD turbulence 

The instability triggered by the super- Alfvenic flow of cos- 
mic rays upstream of a shock has been analyzed in two 
ways: one is related to the reson ant interaction of the cos- 
mic rays with the Alfven waves ( McKenzie fc Vol^ Il982^ 
and is essentially described by a kin etic t h eory. The other 
one has been recently proposed by iBelll l)2004fl and em- 
phasized the importance of non-resonant interactions, in 
which the DC-electric current of cosmic-rays generates a 
Lorentz force responsible for the amplification of the MHD 
perturbations. In fact, this is the return current in the 
background plasma which generates the perturbations un- 
der some conditions. Both resonant and non-resonant in- 
teractions are actually two regimes of the same streaming 
instability. 



2.1. The non-resonant regime of the instability 



The idea developed by IBelll l|2004D states that the 
cosmic-ray fluid weakly responds to perturbations of 
wavelengths shorter than their Larmor radii, so that the 
main response is in the form of a perturbed current of 
the background plasma. Actually, the major role played 
by the cosmic rays is to generate a DC return current in 
the plasma. Because there are resonant interactions with 
cosmic-rays of Larmor radius tl at MHD scales k such 
that kri^{e) — 1, the validity criteria for the dominance of 
non-resonant interactions needs to be analysed carefully. 
Indeed, if one states that it holds for wavelengths shorter 
than the shortest CR-Larmor radius, then one has to pay 
attention to the possibility of going beyond the validity 
of MHD description, requiring fc(5o < 1 where 5q = Va/oJci- 

Let us first reformulate the calculation performed by 

iBelll |2004) as follows. The plasma remains locally neu- 
tral, so that the electric charge carried by the cosmic rays 
is balanced by an electric charge carried by the back- 
ground plasma; no being the number density of electrons 
or protons, the number density of protons in the cos- 
mic ray component is Xp"-o (Xp < 1) and the electrons 
contribute to the CR-population with a number density 
Xeno (xe < 1); the charge density in the CR-population 
is therefore pen = (Xp ~ Xe)'noe. Similarly, the electric 
current generated by the Fermi process upstream is bal- 
anced by an electric current in the thermal plasma; thus 
we state that pcR + Ppi = and Jen + Jp\ = 0, and 
•^CR = iXp ~ Xe)"-oeV^h with respect to the upstream 
medium, Vsh being the shock velocity. Therefore, the elec- 
tron component of the thermal plasma drifts with respect 
to the ion fluid with velocity Vd = ^f^^Vsh- This drift 
does not destabilise the slow magneto-sonic waves as long 
as it remains smaller than the sound velocity Cg, which 
holds when the CR-population is sufficiently tenuous that 
(Xp-Xe)/(1-Xe) < M^^ {M = Vsh/cs is the shock Mach 
number). When the drift velocity exceeds the instability 
threshold, the fast growth of magneto-sonic modes gen- 
erates an anomalous resistivity that diminishes the cur- 
rent down to a value close to the threshold value. Or in a 
more explicit way, the magneto-sonic waves are amplified 
at a rate proportional to Vd — Cg when the drift velocity 
of the thermal electrons exceeds the sound velocity. The 
electrons are scattered by these waves through Landau 
resonant interactions and thus undergo an effective colli- 
sion frequency proportional to the energy of the waves. 
This gives rises to an anomalous resistivity that tends to 
reduce the electric current as the wave amplitude grows 
and thus to make the drift velocity Vd decrease to a value 
close to the sound velocity. 

Hereafter, we assume for simplicity that the magnetic 
field lies along the shock normal. The generalization of our 
results to oblique situations is straightforward, as long as 
a de Hofmann- Teller transformation is possible. 

In the non-resonant regime of the streaming instability, 
the cosmic-ray fluid is, in a first approximation, passive: 
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only the thermal plasma responds while experiencing a 
Lorentz force due to its charge ppi and its current Jpi, 
under the frozen in condition of the magnetic field, namely 
E + u X B = 0; where u is the local fluid velocity. 
The Alfven wave equation is then modified as follows: 



& 



52 



d 



d 



a¥ oz'' po ot oz 



(1) 



This leads to two dispersion relations for right and left 
modes, namely: 



klV? 



(2) 



where 



(3) 



and pq — nip uq. 

The magnetic field Bo and the Alfven velocity Va — 
BoI^Pq Po are to be considered as mean values. The 
scale must be compared to the minimum MHD scale 
So = Va/lOci — c/cjpi, below which MHD no longer applies, 
and one gets fcci5o = \xp — Xe|^sh/VA. From the previous 
dispersion relations, we easily deduce that the waves are 
stable when T4h < Va, and that they become unstable 
only when Vsh is sufficiently larger than Va and fcy < kc- 
Indeed for Va/V^k ^ 1, one of the branch is unstable when 
oj^ — —Vl{k\\kc — /cp < 0. This unstable mode is of right 
or left circular polarization depending on the main compo- 
sition of the CR-fluid and the orientation of the magnetic 
field. Let b be the unitary vector that points toward the 
same direction as the vector —{xp — Xe)Bo] then one gets 
u = ib X u: the mode is thus of right circular polariza- 
tion with respect to the direction defined by b. For the 
likely case of a proton dominated CR-fluid, the mode is 
left-handed with respect to Bo] in other words, it rotates 
in the same sense as the protons. Such a left mode ex- 
ists only for kSo 1, otherwise it is heavily damped by 
resonant cyclotron absorption. 

Modification of the instability due to the CR-response : 
Following Bell's formulation, the modification is described 
by the complex factor cr* (k) such that 



OJ 



k?,Vl±Vlkck\i{l-a^) =0 



(4) 



where corrections in V^/V^ were neglected and cr* reads: 
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(5) 



In this work, we assume a CR distribution function f{p) oc 
p-4-e ]3,gtween po and Pmax, where e may be either pos- 
itive or negative (see paper II, Marcowith et al. 2005); 
moreover 



(6) 



with A = [fc||rL(p)] and A, = (fcy?'*) ^ where 
r* = rL(p»). Here (r*) sets the minimum momentum 



(Larmor radius) of the cosmic-ray distribution function. 
This cut-off depends in principle on the distance x to 
the shock front (measured along the shock normal) , since 
the cosmic ray distribution function roughly decreases 
as f{p,x) cx exp[-a;/£D(p)]; ^d(p) (l/3)crs(p)c/Kih is 
related to the scattering time Ts{p) and is an increasing 
function of p. Hence at each distance x there exists 
p, (x), defined by a; = ^d(p*), such that the contribution 
of the smaller energies p < p* is negligible, since the 
corresponding diffusion lengths are short. The cosmic-ray 
density is tick — J^"""' d^pf{p) is related to the CR 
pressure at the shock front via ncR — 3PcR/i^Pac), with 
$ a dimensionless number of order log(pniax/'7ipc) ^ 10. 

For short waves, k ::§> 1/r* or A, ^ 1, die{(j^} ~ 
and the previous result of Eq. Q holds. In particular, the 
non-resonant growth rate is 

Gn-rcs(fc||)~l^A(fccfc||)'/' (7) 

for l/r» <C fc|| ^ fcg. Note that this non-resonant insta- 
bility is not operative for a wavenumber at distances 
X <C ^d('^l = l/fc||) since the corresponding (x) <C 
r* being an increasing function of x. The exact spatial de- 
pendence of the growth rate will be specified further on. 
The cut-off wavenumber kc is redefined from the Eq.lP)l 
by: 

47rncReysh IStt Pgr ^sh 1 



k, = 



Bo 



^ Bl c 



(8) 



For long waves, fc ^ 1/r, or A* 3> 1, 3fie{(T*} = 1 and 
the CR-response dominates, i.e. the non-resonant insta- 
bility is inactive. 

2.2. The resonant regime of the instability 

The imaginary part of cr» describes the resonant inter- 
action between cosmic rays and Alfven waves; it is re- 
sponsible for a growth rate that reaches a maximum for 
A* = 1, and for longer waves (A, > 1), ^mjcr*} = ||-^t^. 
However we will adopt a slightly different description, in 
the sense that we expect to get oblique Alfven waves that 
essentially are of linear polarisation, which changes the 
resonance conditions as both electrons and ions, moving 
forward or backward, can resonate either with forward 
modes or backward modes. The small instability growth 
rate is the same, within an angular factor of order unity, 
and is given by 

Grcs(fc||,2:) = Go(fc||,x)0(x/^D(fc||)), (9) 

where ^D(fc||) should be understood as ^d(?'l = ^/k\\); the 
exact spatial dependence of (f) will be specified further on. 
The growth rate Go is given by 



Go(fc||,a;) 
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So "•0 V I cos^l ' °" 3 
this expression can be found in iMelrosd l|l98i 
ficient 



Va ] (fc||r*)^- 
the coef- 



ao(£) = i(l + £)(4 + e) 
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Eg. II ion scales as B ^^^^^ similarly to 

IPtuskin fc Zirakshvilil l|20n,'^ . It clearly shows that 
only modes propagating forward are destabilized when 
Vsh is sufficiently larger than the Alfven velocity. The 
resonant growth rate is maximum at the scale r* and 
scales like ri+^ in regards to the distance to shock front 
X. Note that the backward waves arc damped at the same 
rate than the forward waves arc amplifie d. Th e original 
calcul ation has been done by iLerchel ((^^, IWentzei 
lll969tl and use d in the theory of cosmic ray transport by 
ISkilfind (dlzi, then bv lMcKenzie fc Vo lk (1983) for the 
excitation of turbulence upstream of a shock. Hereafter 
we will assume for simplicity e = 0, the value of this 
parameter will be discussed in paper II. 

The function stems from the solution of the evolution 
equation 



5 ^ ^ ^ 

ox ox ox 







(11) 



with (/)(0) = 1. This function (f) represents the spatial 
profile of the CR distribution function, which decreases 
with the characteristic scale given by the diffusion length 
= D/Vsh_- The diffusion length depends on the Larmor 
radius tl, and thus the instability growth rate G'(fc||) de- 
creases with (j) with a characteristic length which is the 
diffusion length for tl — This can be derived rigor- 
ously from the general expression of the growth rate that 
involves a resonance for fc||rL/i = 1 (/i is the particle pitch- 
angle cosine). In the case of uniform diffusion coefficient 
D, (j)^ exp{~x/£B)- 



3. Saturation mechanism and stationary spectra 

3.1. The main elements of the theoretical description 
WKB-approximation : 

The turbulence spectrum is not homogeneous but grows 
when approaching the shock front like the cosmic ray dis- 
tribution. The scale of spatial variation is given at each en- 
ergy by the diffusion length £d = (l/3)(c/V^h)cTs, which 
is much larger than the Larmor radius at the same en- 
ergy; Ts > tL is the scattering time defined further below 
as a function of the turbulence spectrum. Since the mode, 
that undergoes a resonant interaction at this energy, has a 
wavelength equal to 27rrL, its wavelength is shorter than 
the diffusion length and thus a WKB description of the 
evolution of the spectrum is suitable. This statement is 
also true for the non-resonant modes. 

Advection versus non-linear coupling: 

For each spatial scale (or each value of fc||), there are 
three relevant time scales: i) the growth time scales of 
the resonant and non-resonant modes Grcs{k\\T and 
Gn-icsik\\,x)^^ , ii) the non-linear time scale that can be 
defined as the eddy turn over time Tn-\in{k) — (fcu(fc))~^, 
where the turbulent velocity u{k) is such that pQu{k)^ is 



the turbulent energy density at that scale, or by an appro- 
priate non-linear scattering time; iii) the advection time 
Tadv (k) at which the mode is caught up by the shock front 
that propagates faster than the forward waves {Vsh > Va). 



The time scale radv(fc|| 
and one obtains 

1 c2 

Tadv = oT72"'^s(l/^l|) ■ 



£u{k\\)/Vsh with fc||rL 



(12) 



The pitch angle frequency — is known for an 
isotropic spectrum S{k), and, if furthermore the spectrum 
is a power law of the form vka^lni^ / ^rnin)~^ , then 

~ 7r{l3 - 1)lo^vP^-' , (13) 
with the rigidity defined as p = femin'^L [see ICasse et alJ 



l|2002|) ]. In the above equation, the prefactor (3—1 should 

be replaced by [log(fcmax/fcmin)] if /3 = 1. The turbu- 
lence level rj is defined in the next paragraph. 

However, we deal with anisotropic spectra of the form 
5'3d c>c kj^k^'^ with q > 2, leading to the same formula 
Eq. H13|) (see paper II, for further details). Hereafter, we 
will use ID-spectra S{k^^) defined such that 



which implies notably 



5(fc||) , 



(14) 



(15) 



For convenience, hereafter, the normalization of S is de- 
fined in regards to the magnetic energy density at the 
infinity 



dfc|| 
"2^ 



S{k\ 



dlog(fc| 



SB^{k«) 



Bl 



SB^ 



(16) 



where SB is the turbulent field amplitude, and Boa the 
(original) uniform component to be taken in the interstel- 
lar medium far upstream of the forward SN shock. The 
Alfven velocity Vaoo in the interstellar medium can then 
be deduced immediately. The magnetic field and Alfven 
velocity amplified by the streaming instability are here- 
after noted B and Va- The two Alfven velocities are 
linked by Fa — V\na/(l — vV^'^ [see lPtuskin fc Zirakshvilil 
(Hoal]- The quantity t] = 6 B^ / {S B"^ + B"^) determines the 
strength of turbulence; in particular 77 ^ 1 corresponds to 
5B/ Boo — > -|-oo. Finally, the magnetic field turbulent am- 
plitude at a scale fci| and the ID spectrum S'(fcii) are tied 
by the relation: 5B^{k{) = fc||S'(fc||)/(27r). 

While considering the resonant instability, we distin- 
guish the forward ID-spectrum of forward waves and 
the spectrum S~ of backward waves. The pitch angle fre- 
quency is the sum of the two contributions because par- 
ticles resonantly interact with both spectra irrespectively 
of direction of motion (this is an important point related 
to the resonance condition with linearly polarized Alfven 
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waves, as mentioned before). The advection time is using 
the definition of in Ea. H13() 



T"adv(fc|| 



5- 1 - 77 ■ 



(17) 



Note that the advection time remains finite when the mag- 
netic field tends to be completely turbulent, i.e. when 

We introduce an important dimensionless quantity 
that measures the ratio of advection time to instability 
growth time : a(fc||) = 2Gradv For the resonant regime of 
the instability using Eq. ifTUI) and Eq. lfT7|l we obtain: 



7r X B 



(18) 



where Maoo ~ V^h/VAoo is the Alfvenic Mach number 
measured with respect to the interstellar magnetic field 
value, ^cR = Pcr/poV^ is the ratio of CR to shock pres- 
sure and can re ach values as high as . 5 in non-linear accel- 
eration models ijBerezhko et a l ll996HBerezhko fc Elhsonl 

Eiai), and b'^ = SB^ + B^. 

The ratio B/B^o stems from the spatial dependence 
of in Eq. (|10|) and from the 1 — 77 factor in Eq. H17|) : 
B/Boo depends only on the distance to the shock front. 
The Eq. H18|l accounts for the amplification of the mag- 
netic field along the normal to the shock front and permits 
the inclusion of both non-resonant and resonant regimes 
in the evolution equation (see Appendix A and B). We 
define the reference spectrum 5'*(fc||) 



'k\\ Br, 



(19) 



Simultaneously the Alfvenic turbulence develops an 
energy transfer mainly in the transverse direction which 
determines the shape of the transverse spectrum in k±. 
The non-linear transfer rate is t~}_ii^ «:! fc_L m_l ~ Vaoo x 

^^5*30(^11 , fcj_)]^^^. We can define the efficiency of the 
energy transfer process using Eq. IjlGI) and the dimension- 
less number ka 



KA 



'^adv 
'^n— lin 



cVa 1 - ?? 



(20) 



This number is sufficiently high for the Alfvenic cas- 
cade to fully develop; this will be discussed in Section rOl 
Because the Alfvenic cascade does not convert energy from 
the forward waves into backward waves, the backscatter- 
ing of Alfven waves off slow magneto-sonic waves will also 
be considered and will proved to be efficient to redistribute 
the energy over all the spectra. This discussion is post- 
poned to Section^ 

3.2. The spatial profiles and spectra 

As mentioned earlier, the quantity r^{x) denotes the mini- 
mum Larmor radius of streaming cosmic rays at a distance 




Fig. 1. Sketch of the evolution of the turbulence spectrum 
5' in the (x, A:||) plane (upper graph) and of the cosmic-ray 
distribution function f{x,p) as a function of x for two 
values of the momenta (lower graph). See text for details. 



X from the shock front; 
dition: 



da;' 



fe[a;',r*(a;')] 



= 1, 



,{x) can be defined by the con- 



(21) 



which, if Id does not depend on x, amounts to a; = 
iu [r*{x)]. 

The non resonant regime of instability occurs for 
modes such that fc||r, (a;) 3> 1, hence at distances x > 
a;min(fc||) with Xmin(fc||) defined by r* [a;,nin (fc|| )] = l/fc||. Of 
course, this non-resonant growth occurs provided there ex- 
ists cosmic rays with tl > r*{x), hence for x < x^ax ~ 
^D(''L,max)- In Contrast, the resonant interaction growth 
rate is maximal for fc|| ~ 1/r, (x), therefore the vicinity of 
the shock front x < x^-nin{k\\) is dominated by the resonant 
regime of the instability. 

In Fig. we sketch the simultaneous evolution of the 
turbulence spectra as determined as a function of x and fc|| 
and the distribution function of cosmic rays as a function 
X for two values of the momenta (hence two values of the 
Larmor radius). At point X2, far from the shock front lo- 
cated at a; = 0, only the non-resonant instability has been 
active in the momentum range r~^{x2) < fc|| < kc{x2) 
and has produced a spectrum S cx fc|^^ (Section 13.2.1(1 . 
At point X2 as there are no cosmic rays with Larmor 
radii such that the resonance condition k^^ri^ ^ 1 can 
be satisfied in the above momentum range, since r^(x2) 
is by definition the minimum Larmor radius of cosmic 
rays at point X2. The lower graph of Fig. sketches 
accordingly the evolution of the distribution function of 
cosmic-rays at momentum Pi,(x2), which corresponds to 
the minimum momentum of cosmic-rays at X2, or, equiva- 
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lently, to a Larmor radius rj,(x2). The minimum momen- 
tum (or Larmor radius r*) is a growing function of 
in particular Pi,(x2) > J5*(xi). Hence, at point xi, the res- 
onant instability has been active in the momentum range 
fc|| < l/r^(2;i) as there are cosmic-rays with Larmor radii 
that can satisfy the resonance condition in this range. The 
corresponding spectrum is S* oc k^^ (Section I3.2.2() . At 
momenta k\\ ^ l/r-i,{xi), only the non-resonant instabil- 
ity has been active, for reasons similar to those discussed 
above, and therefore the spectrum S oc Finally, at 
the shock front, the resonant instability has overtaken the 
non-resonant instability over all the wavelength range, so 
that the final spectrum S oc k'^^ . 

Since the two regimes of the instability take place in 
different regions, we can solve for the spectra using in 
a first place the equation involving solely non-resonant 
growth for x-cain{k\\ ) < x < Xmax and then use this solution 
at a;min(fc||) as the initial condition for resonant growth up 
to the shock front at a; = 0. 



saturation mechanism occurred earlier. Even if we ac- 
counted for a variation of 77, the amplification level would 
still be too large, as the following saturation mechanism 
keeps the magnetic field energy density to a lower level. 

Saturation mechanism : 

In fact, the non- resonant growth should saturate 
much earlier. As the magnetic field gets amplified 
beyond its initial value, one can extrapolate the pre- 
vious calculations by substituting B^k for Boo, where 
Bl:k = Bio + /fel„ dlogfc[| 5B2(fc[|) represents the average 
field on scales larger than fc^^, and again 5B{k\\) denotes 
the amplified random component on scale fc|| . 

For |xp~Xe| ~ Xpi using Eq. (PJ, the cut-off wavenum- 
ber of the instability kc can be written as 



o„ V \ l/(2-/3<) 



(^min^ 



1/(2-^ /3<) 



(23) 



where /?< is the index of the turbulent spectrum in the 
vicinity of the shock front, i.e. that which is produced by 
the resonant instability. Introducing 



^ ~ (37700)1-"^°^ I c 



12ttP( 



CR 



$ Bl. 



1/2 



(24) 



with m = (3 - 2/3<)/(4 - 2/?<), then 



ATrn^eVsh 127r PcR Vsh 1 



B. 



<k 



$ B%^ c 



(27) 



3.2.1. The non-resonant regime 

The equation governing the growth of the turbulent spec- 
trum through the non-resonant instability is, using Eq.l(7I): 

T/sh^^(fc||,a;) = -2VA^^Jk,{x)^0{xm.^-x)S{k\\,x) , (22) 

where kc{x) is defined in Eq. I|27(l . Assuming for the mo- 
ment that 77 is constant {— ?7oo), Eq. H21|l then gives 
X = ^d(?'*), or: 



The non-resonant instability occurs for wavenumbers 
l/r* < /c[| < kc and its saturation is achieved once 
fccT* = 1. This simple condition leads to a magnetic field 
energy density: 



52 



2$ c 



Pgr 



(28) 



This last estimate is in agreement with Bell (2004). 

However, the instability may also saturate through 
non-linear transfer effects. As the field builds up through 
the instability, the non-linear transfer time tn-iin(fc||) 
along the fc|| direction decreases to the point where the in- 
stability saturates when G'n_res(fc||)in-iin(fc||) = 1. In order 
to see when this happens, one can express the non-linear 
transfer time as: 



<n-lin(fc||) = [k\\VA{k\\)\ 



(29) 



with VA{k\\) ~ 5i3(fc|| )/^/47mo, and the non-resonant 
growth rate G'n-rcs(fc||) = \Jk\\kc V k- We then find for 
the saturated field at a scale k\\ : 



S{kii,x) = S'(fc||,a;i„ax)expJ - — W— L X 

TO V fcmin 



[(^min-^) (^miii'^niax) 

The amplification factor is thus: 

^(fcji , a; = 0) 



127r 



exp 



c 



(25) 



(26) 



It turns out that fcmina;max = ^D(7'L,max)/?'L,max > 1, and 
for /3< = 1 (corresponding to Bohm scaling, see next sec- 
tion) C ~ 1. The amplification would be enormous, the 
instability would deplete the shock quickly, unless another 



<5B^(fc||)^— Pgr 



c kur^. 



(30) 



Integrating this result over fc||r* > 1, we exactly obtain 
the same saturation level as previously, which is quite 
remarkable. This means that even if there is a kind of 
quasi-linear saturation at work, within the same time the 
non-linear transfer remodels the spectrum. The spectrum 
derived above from this remodeling process is likely cor- 
rect but would require a more elaborated theory together 
with sophisticated numerical simulations to be confirmed. 

Using Ea. (|30|l the non- resonant spectrum profile is 
thus 



S{k\i,x) 



27r 6B'^{k\\ 
fc|| 



3^, Vsi^poVl 



1 



pj « -(31) 
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We obtained this result by assuming energy transfer 
along the direction. However the relevant transfer could 
be transverse through an Alfven cascade, as is often con- 
sidered. In this case the transfer rate k±u is supposed to 
be faster and a balance k±u ^ Gn~rcs{k\\) is expected. 
However it generally requires the prescription of a re- 
lation between k± and fc|| [see subsection (|3.3|) ]. Setting 
k±ri, ^ (fc||ri,)™ with to > 1, we obtain a spectrum 



^(fcii) cx kcrl{ki\r^,) 



-2m 



(32) 



There is no clear constraint on the value of to, which may 
be taken to 1 without apparent inconsistency. In any case, 
as to > 1, the same estimate (|28|l of the saturation level 
is found, which makes this estimate fairly robust. 

To draw a complete picture of the saturation spectrum, 
one would need to follow the evolution of non-linear tur- 
bulence transfer at the same time as the evolution of the 
non-resonant instability and the position dependence of 
r*. 

3.2.2. The resonant regime 

As the turbulence is advected to the shock front, cosmic 
rays with Larmor radius tl = l/fc|| appear and induce the 
resonant instability. This latter is then quenched by ad- 
vection through the shock, which thus provides the main 
saturation mechanism of this instability. 

i) For the sake of clarity, we first assume equally am- 
plified forward and backward spectra, and S — + S" . 
The evolution equation for this latter reads 



-x/fd 



-5(fc|| 



(33) 



The initial condition for integration lies &i x = Xmin(fc||) 
and is given by the spectrum produced by the non- 
resonant instability. We argue that the resonant amplifica- 
tion will not modify the overall magnetic field strength B 
over that produced the non-resonant instability by a large 
factor for typical SNe environment and shock wave val- 
ues. Hence, to solve Eq. (|33|l above, we first assume that 
the ratio (1 — 77) ~^ is constant between the shock front 
{x = 0) and the point where the resonant instability first 
comes into play, defined hy x — 2;min(fc||) for fc|| = fcmin- 
Then 1/(1 - r;)!/^ = B^/Boo with B^ the value of B at 
this latter point. In this case, does not depend on x in 
the interval in which the resonant instability acts, and the 
equation for S can be solved as: 



S{kl\,x) = 5(^11, Xniin) + 



-x/tu 



1^ , (34) 



with fires = (7r/<I>)A^AooCcR (see Eq.l^Hll)- The term 
5nr = ^(fcii , Xniin) is the spcctrum produced by the non- 
resonant instability and is proportional to k^'^ due to sat- 
uration. The second term on the rhs is the contribution 
of the resonant instability and we denote it 5*1(^:11) at the 
shock front; it is proportional to fc,?^ and dominates the 



former over the whole range of wavenumbers if aj-csB^,/ Boc 
is sufficiently large. In case of the absence of non-resonant 
instability, = B^o and the expression of arsa lea ds to a 
resonant saturation level fsee lBell fc Lucel3 ()200l|) ') 



B' 
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^crMac 



(35) 



The contribution of the resonant instability to the mag- 
netic field energy density B'^^^/8tt is obtained from: 

dfc|| 5*1(^11) arcs ^ log (^^^) . (36) 



'fcmin 

Hence 

res 

B^_ 



^ 1 



kyi 



Bqc 

Bf: 



sh 



,(37) 



this ratio is larger than one by a factor of a few for shock 
velocities lower than a few times ~ fcRC- Among other 
things, this implies that the scaling with fc|| of the total 
spectrum at the shock front will be dominated by the res- 
onant contribution. 

One can investigate the effect of the above assumption 
1/(1 — 7]) = constant by integrating formally Eq. (|33|l as: 



S{k\i,x) 



^(fc„,..)+^rd.e-^[^, 

k\\ Jo Boo 



(38) 



with y = x/i(i{x)- The integration in the rhs can be un- 
derstood as a function of fc|| so that the dependence of 
the second term on the rhs is not strictly speaking oc fc|^^. 

However B{x)/Boo is a decreasing function of x and there- 
fore it is bounded by below by B[a;*(fc|| = fcmin)]/-Soo and 
by above by B{0)/Bao- The integral in Eq. Ip^ is thus 
bounded by two constants that are independent of fc|| and 
whose ratio is a factor B{0)/B[x^{k\\ = kjnin)] ~ a few. 
Hence the integral modulates only weakly the powerlaw 
cx and we conclude that, at the shock front, the spec- 
trum S cx k'^^ . It should be noted that the total variation 
of the magnetic field in the resonant region is even more 
limited if e > as we shall see in paper II, as in that case 
flres varies as {B/B 

ii) If the backward spectrum is not amplified nor remod- 
eled by mode coupling, the solution is similar to the pre- 
vious case without the factor 2. 

iii) If the backward waves are damped at the same rate 
as the forward waves are amplified, S^S^ = constant. If 
the forward spectrum is amplified an obvious solution is 
obtain for 5"+ S'+(x = a;inin) 

5*+ ~ ^^(fcii) X (exp(-x/4) - e~^) 



s 



s+ 



(39) 



The backward waves are damped exponentially, 
iv) The backward waves can also be generated by 
backscattering process of forward Alfven waves off acous- 
tic waves (slow magneto-sonic modes precisely) . This pro- 
cess deserves a specific development presented in section 
4. 
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3.3. The behavior of Alfvenic turbulence 

The behavior of moderate MHD turbulence - moderate 
in the sense that a significant mean magnetic field is pre- 
served - is peculiar when incompressibility is assumed, 
because of the particularity of Alfven waves dynamics. 
Resonant three wave interactions do not develop as usual 
dispersive waves, because of their specific dispersion rela- 
tion: uj = fc||VA- However, a forward wave and an oppo- 
site backward wave can couple throu gh a resonant interac- 
tion with a third w"a.ve wi th fc|| = l|Bhattachariee fc NeL 
l200lt iGaltier et all liooofl . The weak turbulence descrip- 
tion shows that the energy cascade in the inertial range 
occurs only in the transverse direction to the mean field. 
The stationary spectrum is in kj^ and the dependence 
in fc|| is arbitrary, which means that it is determined by 
the mechanism of generation of the turbulence. This be- 
havior has been observed in numeric al simulations even in 
the re gime of moderate turbulence I Bhattachariee fc Nd 
The extension of the resonant three wave in- 
teraction by taking account of a nonlinear broadening 
due to the relaxation of triple correlation - the so-called 
Eddy Damping Quasi Normal Markovian description - 
ha s been done for MHD turb ulence with a mean field 
by iGoldreich fc SridhaJ l(l99,'il) . They have argued that 
some re-organization of the spectrum occurs in fc|| due 
to scaling constraints between the nonlinear transfer in 
the transverse direction and the parallel propagation of 
Alfven waves. Let us summarize this discussion. These 
Alfven waves (also called shear Alfven waves, as opposed 
to MHD waves that have parallel components) are incom- 
pressible and purely transverse to the mean field. The tur- 
bulent energy density e oc and the eddy turn over time 
■''n-iin ^ {k±u±)~^ . The scale invariant spectrum of the 
energy cascade is unavoidably anisotropic 530(^.1, fc|| ) oc 
k±'^k'^^ . The critical balance assumption of Goldreich and 
Shridar is that the transfer rate Ttr ~ Tn_iin ~ ta at all 
scales, where ta = (fc|| Va)"^. Then for an anisotropic iner- 
tial cascade such that the energy transfer rate at each scale 
is constant, namely Q oc e/Tn-vm oc k±u^j_ — constant, a 
relation between parallel and transverse wavenumbers is 
found: 



fell 



)l/3 



Va 



2/3 



(40) 



In this more elaborated description (EDQNM) no en- 
ergy transfer from forward waves to backward waves 
and vice versa takes place. Recent numerical simulations 
l)Cho fc Vishnia c. 2000; Maron fc Goldreich, 2001} have 
suggested that the scaling T^-nn oc ta was preserved in all 
regimes, so that Tn-un = ta/x, with x a constant indepen- 
dent of the wavenumber at all scale, and Ttr = Tn-Vm/x- 
The previous relation is thus extended to 



Q 



1/3 



i.2/3 



(41) 



Whereas the weak turbulence theory leads to a spec- 
trum in A:J^/(fc||) with an arbitrary function / of fcy. 



the Goldreich -Shridar theory leads to a spectrum S^b oc 
^-'3-2/3j^^^^^^2/3^^ When a scale invariance in fc|| is gen- 
erated in the turbulence situation, the spectrum is of the 
form 5*30 oc kj^'^k^^ . Then Eq. together with the as- 
sumption of a constant energy transfer rate Q oc e/Tn-Vm 
with e ~ k'^k^^S^Y) provide a relation between the index of 
the parallel spectrum with the index a of the perpendic- 
ular spectrum (Ga ltier et al... ,200.5.) : 



3a + 2/3 = 7 with a = q-l . 



(42) 



This is considered t o be the generaliz a tion of Iros h nikov- 
Kraichnan theory l)lroshnikovl Il964t iKraichnanl Il965|) 

when anisotropic effects are taken into account. Still some 
arbitrariness is maintained. However the CR-instability 
in resonant regime generates a turbulent spectrum such 
that /3 = 1, and the transverse Alfvenic couplings be- 
tween modes then lead to a = 5/3. Only couplings with 
slow magneto-sonic modes may allow to obtain the same 
spectrum for the backward waves and the slow waves. 

4. Nonlinear generation of backward waves 

The process yl+ ^ + 5+, where A represents Alfvenic 
modes and S a slow magneto-sonic mode, is the only 
process that can transfer energy from forward waves to 
backward waves and it turns out to be efficient, as will 
be seen further on. The frequency of the slow magneto- 
sonic mode is such that lu" — fc'*V^,n(6's); for convenience, 

we write it = /3sfc|jVA, where, for Cg < Va, /3s — 

2 

■^(1+-^ sin6's)~^"; this number is assumed smaller than 
unity and weakly varying with Og ■ The process is most ef- 
ficient under the resonance condition : uj^ — lu^ — uj" = . 
Since the wave vectors are such that — k — = 0, 
we obtain the following relations between the parallel 
wavenumbers: 



'II 



1- A 



< . 



(43) 
(44) 



Therefore, when the magnetic field is above the equipar- 
tition value (e.g. Va > Cg), /3s < 1, and we always 
get a backscattering of Alfven waves off slow magneto- 
sonic modes. Backscattering would not be possible with 
other MHD waves, for obvious kinematic reasons. This 
backscattering process with Alfven waves is analogous to 
the Brillouin backscattering process with usual electro- 
magnetic waves of the vacuum. Even if no sonic waves 
are excited beforehand, the primary Alfven waves can 
gene rate them spo ntaneously above some threshold [see 
■Pelletier fc Kersala (2000)]. In the interstellar medium 
Va — 3 Cg; it is already sufficient to get the backscattering 
process. The domination of the Alfven velocity over the 
sound speed is even increased at the external shock of SNr 
because of their convexity. Indeed in the shock frame, the 
ambient medium converges towards the front at a veloc- 
ity — Vsh that points towards the center of curvature, the 
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density increases and therefore the frozen in magnetic field 
has an amplified transverse component. From the evolu- 
tion equations multiplied by T^dv: we introduce the dimen- 
sionless parameter 



TT cVa 



(45) 



which measures the importance of the backscattering pro- 
cess as compared to advection. Typically, k is a number 
close to one. 

For the most interesting cases where the backscattering 
process efficiently remodels the stable spectra, the asymp- 
totic spectra, determined externally by the turbulence 
in the interstellar medium, can be ignored. The spec- 
tra are then proportional to S'*(fc||)A:]^'. Because the dif- 
fusion length is generally not constant, but dominated 
by the spectrum of unstable waves, for numerical simu- 
lation purpose, it is convenient to describe the profiles 
of the wave spectra with the help of a dimensionless 
variable y defined by dx = £d(^l = ^/k\\,x)dy. Then 
the function (f> = e^^. We have to bear in mind that, 
when the problem is solved for the variable y, we can 
reconstruct the spectrum profile in the variable x. Since 



iuirL - l/fc|| 
one finds: 

y = • 



,X) - {l/3)ic/Vsh)k7\kJknuny-^V'Hx), 



.Vsl 



2-/3 



The quantity r, (x) is defined hy y — 1 for tl — r^,, hence 
y = (fciir,)^ . The regions with y > 1 {y < I) correspond 
to far (close) distances to the shock front and is dominated 
by the non-resonant (resonant) waves. Therefore the evo- 
lution of the spectra reduces to a differential system that 
governs the evolution of their amplitude as a function of 
the y variable. As long as k is small, the solution given 
by Ea. H39|) is slightly modified and the order three sys- 
tem that describes the generation of backward A-wave 
and forward S-wave is sufficient (see appendix A). The 
numerical integration leads to the solutions sketched on 
figure 121 that shows that for, increasing k, more and more 
conversion into backward A-waves and forward S-waves 
is realized. However when k is increased significantly, one 
has to take account of the secondary process where back- 
ward A-waves decay into forward A-waves and backward 
S-waves. The evolution is then described by a system of 
order four (see appendix B). The numerical solutions are 
displayed on figure|3| It can be that a significant backward 
spectrum is generated; however without changing the or- 
der of magnitude of the primary spectrum. 

5. Downstream: dynamo action and turbulence 
relaxation 

If turbulence is still moderate downstream, then the spec- 
tra built upstream are transmitted across the front, and 
thus a k'^^ ID-spectrum is maintained downstream. Bohm 
diffusion would then applies downstream as well. 



X 



5 7" 
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A* 
A" 
S 



>- 
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_o 

-10 



K = 0.5 

X„ = (BoVBj)- 
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A" 
S 



2 4 6 8 10 

y 

Fig. 2. Solutions of third order system for different value 
of the K parameter (see the definition in the text). Even 
if they have been defined for y between and 1, the reso- 
nant waves profiles are calculated between the shock front 
and y = 10 where they match their asymptotic interstellar 
values. The wave-particle resonance depends on the par- 
ticle pitch-angle a, at tl = we have k\\r^cosa ~ 1, 
the product k\\r^ and then y can be above 1. In the up- 
per panel, k = 1 is high leading to a strong conversion 
of forward Alfven waves into backward Alfven waves and 
sound waves. For y < 5, the resonant instability takes over 
the non-linear transfer, the forward Alfven waves are pro- 
duced and the backward waves are pumped. The sound 
waves are heavily produced between y = 5 and the shock 
front. The ratio of forward to backward Alfven waves at 
the shock front is about three orders of magnitude. In the 
lower panel n = 0.1, the production of backward Alfven 
and sound waves is less intense. In both cases the amplifi- 
cation factor B:,/Boo — 10. The boundary conditions are: 
X+ = X- = Xoo = A-"^ < 1 and X' = 0. All simulations 
have been performed with e = 0. 



5.1. Helicity and estimate of dynamo amplification 
downstream 

The non-resonant regime of the streaming instability in- 
duces a left-right symmetry breaking. Therefore the tur- 
bulence carries helicity which offers grounds for dynamo 
action. The helicity can be calculated in term of the dif- 
ference between the spectrum of right-handed modes Skh 
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Fig. 3. Solutions of the fourth order system. The non- 
linear transfer parameter n = 1 and the damping rate of 
sound waves is = 5 x 10""^. This choice corresponds to 
a ratio Va/cs = 3. All simulations have been performed 
with e = 0. The results show that as long as the damp- 
ing or growing effects do not dominate the transfer among 
the waves (for y < 5), the non- linear transfer mostly pro- 
duce backward and forward sound waves equally. In both 
cases the amplification factor B^,/Boc = 10. The bound- 
ary conditions are: = = Xoo = A^"^ <C 1, and 
=x'- = 0. 

and the spectrum of left-handed modes S'lh (the k± de- 
pendence has been integrated out): 

(5rh - 5LH)fc|| . (46) 

The integrand can be considered as the helicity spec- 
trum 5*^- This spectrum is used to calculate the so-called 
"alpha" -parameter of the turbulent dynamo theory: 



ttD 



r(fc|| 



r2(fc|| 



:'5'if(fc||)' 



27r 



(47) 



where F is the damping rate of the turbulence in stationary 
state. In our problem the main damping mechanism is the 
shock advection: r(A:||) = 1/Tadv(fc||)- For the non-resonant 
modes ajf <C F2(/c||) and the dynamo coefficient reads: 



2 VI 
3 



5, 



RH 



'LH 



1 dfcii 



fcn 2tt 



2c Vl , r 
4-ln- 



ro 



(48) 



The helicity is t ransfered th r ough the shock as has 
been calculated by ISchlickeise^ l)l998(l . Helicity in the 
spectrum matrix leads to a third diffusion coefficient for 
the cosmic rays because the two transverse space variable 
are correlated (Aa;iAa;2) ^ (see paper II). 

The mean field evolves in the turbulent plasma accord- 
ing to the following equation: 



d 

— A = a-oB 
at 



- u X B + j^t AA 



(49) 



where i>t is the turbulent magnetic diffusivity. A typical 
scale for the variation of the mean field arises, namely 



idyn — vt/a^ with an associated time scale Tdyn = 
f-dyn/u2- Morc precisely, the dynamo modes of wavelength 
larger than Vt/otTi grow and it is expected that the mean 
field reaches an intensity on the order of the equipartition 
value, not more. 

5.2. Relaxation or compression downstream 

The turbulence properties downstream (level of turbu- 
lence, spectral index) can be constrained from the size 
of the X-ray filaments in young SNr [see IParizot et al. I 
(2006)]. It is shown that the relativistic electrons with tens 
of TeV energies producing the observed synchrotron radi- 
ation hav e a diffusion coeffi cient close to the Bohm value. 
Recently, IPohl et all (jiflQ^ have pointed out the impor- 
tance of the relaxation length of the turbulence down- 
stream the shock. The authors stressed that the size of 
the X-ray filaments observed in young SNr (see the dis- 
cussion in section may be well controled by the tur- 
bulence rather than by the synchrotron losses. However, 
the previous analysis has been made assuming an isotropic 
turbulence spectrum, which is not correct at least for two 
reasons: the turbulence is already anisotropic upstream 
and the magnetic field amplified upstream is compressed 
in the direction parallel to the shock front. 

In order to elucidate the way the turbulence acts on 
relativistic particles, we compare the non-linear Alfven 
transfer time tn-iin(fc||) = [^||^A(fc||)]~"'^ 

t„-Un(fc||) = - ^ (/||fc||)(/^-3)/2 , (50) 

and the term /3 — 1 in the prefactor should be replaced by 
1/ log(fcniax/fcmin) whcu (3 = 1. Thc dowustrcam return 
timescale reads: 

tret ^ TT ] ots, (51) 



ts = {2n)~^?]~^{l\^/c){ri^/l^^)^~^ being the angular scat- 
tering timescale downstream; in the above expressions, 
Z|| is the coherence length in the parallel direction. The 
prefactor k < 1 accounts for the shortening of the re- 
turn timescale in compressed turbulence. The condition 
^rct < in-iin mcaus the particle does explore distances 
smaller than the relaxation length of the turbulence down- 
stream and experience a compressed rather than a relaxed 
turbulence during their journey downstream. In the case 
of Bohm type turbulence (meaning /3 = 1) the ratio of the 
two timescales translates into 

^ret ^ K C log(fcniax/femin)~^^^ 



(52) 



Using typical values of the magnetic field, mean 
density and shock velocity in our problem (see next 
section) and acknowledging for a saturation level of 
Maoo(,ckB*/Boo ^ ■n/{l ~ v) - 1/(1 - V), we find a 
ratio irot/in-iin ^ 0.5k (V;^;LiCcr/'J')^/'* where V^h -i is 
thc shock velocity in units of 0.1c. The previous ratio is 
thus expected to be < 1 unless the shock velocity is lower 
than lO^^c . 
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6. Astrophysical consequences 

Supernova blast waves explore different external interstel- 
lar media (ISM) during their evolution. We consider only 
two phases of the SN evolution: the very early free ex- 
pansion phase where the shock velocity can reach extreme 
values as high as T4h — 0.1c and the late free expansion 
phase (or early Sedov self-similar phase) where the shock 
velocity drops to 14 ~ 10~^c. These two phases are the 
more relevant concerning high energy cosmic ray produc- 
tion. In the latter phase the remnant may expand either in 



a hot rarefied (T ~ 10^ K, 



10 ^ cm ^) interior of 



a massive star wind bubble or in a warm partially ionised 
(T ~ 10'' K, n ~ 10-1/° cm-3). The mean ISM magnetic 
field in both cases is conservatively taken to 3/xGauss. In 
the very early phase of the SNr evolution, the medium is 
probably much denser with n ~ 10 — 100 cm"'^. 

The ratio of the saturation magnetic field energy den- 
sity of the non-resonant and resonant instability (when 
both regimes are present) respectively is S^^ycs/ S^cs — 
^/VsiJlXcB~cj ■ The non-resonant regime appears to dom- 
inate for the ve ry ea r ly fre e-expansion phase as already 
pointed out by iBelll l)200-d) . while the resonant regime 
dominates by a factor 10 — 100 (see discussion in sec- 
tion I3.2.2|) in the late free-expansion phase and Sedov 
self-similar phases. The magnetic field deduced (only lower 
limits) from the siz e of bright X-ra y filament s in young SNr 
llBerezhko fc VolkL l2004HVink L |2004; Vo lk et all 120051 
|Parizor'et'^n ~ l2006(l is expected to be mostly produced 
in the resonant regime should then scale approximately 
as VVsh [see Eq.(jSHl)]- We have seen above that if the 
non-resonant regime contributes substantially to the am- 
plification this dependence is not as simple and the way 
the shock decelerates during the earlier phases can modify 
it. The amplification by the non-resonant instability may 
lead in the most extreme cases to very high amplification 
levels, pointing towards SNr in very early free-expansion 
phase as efficient CR accelerators. In that case, the i nag- 
netic density should scale as as pointed out bv iBelll 
l)2004|) . This issue is of prime importance and should de- 
serve detailed observational investigations, unfortunately 
difficult to perform in this SNr evolution stage. However, 
this early phase lasts for a very small fraction of the whole 
SNr lifetime except in a low density and highly magne- 
tised medium as expected in a turbulent hot ISM phase of- 
ten ca lled as superbubbles [see for instance IParizot et al. I 
llo^)] . Answering the question of the maximum CR en- 
ergy expected in SNr and the origin of the CR knee at 
'--^ 3 X 10^^ eV requi res then a time dependent C R spec- 
trum calculation fsee lPtuskin fc Zirakshvilil l| 2005')] and to 
account for the CR diffusion regimes correctly. If the first 
point is beyond the scope of the present work, the second 
point will be discussed in paper II. 



7. Conclusion 

We summarize the main results of the work as fol- 
lows. Upstream of an astrophysical shock, the cosmic ray 



streaming triggers an instability that has two different 
regimes: one occurs under resonant condition and dom- 
inates the longer wavelengths of the Alfven spectrum, the 
other occurs off resonance and dominates at shorter wave- 
lengt hs [this is the Bell regime of the instability ijBellL 
In the purpose of investigating the turbulent trans- 
port of the highest energy cosmic rays both regimes have 
to be considered over the remnant evolution. The non- 
resonant instability saturates either by non-linear transfer 
effects or by a quenching effect at kcr^. = 1. The satura- 
tion level is fc||5(fc„) ~ ^cR{Vsh/c) M^^/($/c||r,). The 
main saturation mechanism for the resonant instability 
stems from the fact that the shock front catches up with 
the growing waves over a diffusion length. The saturation 
level is modulated by the non-resonant saturation level 
according to: A:||5(/e||) ~ ^CR (A^Aoo/*) x Bn-ros/^oo- 
The non-resonant regime of the instability dominates the 
resonant contribution only for very fast shock velocity. 

The streaming instability partially determines the 
spectrum, namely its fcy dependence. The spectra are 
close to which, as will be shown in the second 

paper, can lead to a Bohm scaling for the transport 
of cosmic rays. The k± dependence of the spectrum is 
remodeled by the non- linear cascade of Alfven waves, 
that essentially works transversally, the transfer time 
being short enough as compared to the advection time. 
The excited Alfven turbulence constitutes the scattering 
medium for the cosmic rays and it would be incomplete 
if only the forward Alfven waves would be present as 
a result of the resonant streaming instability. A second 
nonlinear transfer develops which is the backscattering of 
primary Alfven waves off slow magneto-sonic modes, and 
which re-distributes the energy from the forward Alfven 
waves to the backward ones and to the magneto-sonic 
ones. This process turns out to be unavoidable because 
the Alfven speed exceed the sound speed upstream and is 
sufficiently fast compared to the advection time. 
Some turbulent dynamo action can be expected down- 
stream, but the intensity of the mean field should not 
significantly exceed the equipartition value. The turbu- 
lence is also compressed at the shock front producing 
reduced residence time of the relativistic particles down- 
stream. Apart for low shock velocities, i.e. lower than 
IQ-^ c, the residence time downstream is lower than the 
non-linear transfer time controled by the Alfvenic cascade. 
The turbulence spectra downstream are likely similar to 
those that have been formed upstream, since they corre- 
spond to stationary solutions of the turbulence equations 
in a mean field, according to recent developments of 
the theory of Alfvenic turbulence. This implies that the 
Bohm regime of cosmic ray transport if true upstream 
would also apply downstream. Regarding the magnetic 
turbulence in the Galaxy where a mean field imposes its 
constraint, a question rises about the determination of 
the two indices a and (3 of the anisotropic s pectrum. They 
are o nly linked by the relation 3a -I- 2/3 = 7 ( Caltier et all 



l200fil) . It is reasonable to think that /3 and thus a are 
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determined by shocks and thus /3 = 1 and a — 5/3 
(anisotropic Kolmogorov spectrum) would be ubiquitous. 
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Appendix A: The nonlinear operator of 
backscattering 

The energy spectra are normalized such that cj°A^° = 
5*3^7/90 V^; in other words, the N'^ have the dimension of 
an action (namely, an occupation number times H). 



N+ = - 



(27r)3 (27r)3 
X {N+N- + N+N' - N-N') 



, d^k+ d^k" 



(27r)3 (27r)3 

X {N+N- + N+N" - N-N"") 

d3fc+ d3fc- 



(27t)3 (27r)3 
X {N+N" + N+N'' - N-N') 



(A.l) 



(A.2) 



(A.3) 



The transition probabiUt y has been calculated by 
lAkhiezer fc Akhiezed l|l975l) and reads 



w = -/— — — X 

8 po uj+ujfUJs 

(27r)3<5(fc+ -k ~ k')6{tJ+ -uj- - Lj") 



(A.4) 



with the angular factor / depends on the unitary vectors 
n = k/k and defined by 



(n|)2(nf)2 



(A.5) 



where ei, is the unitary vector in the direction of the mean 
field Bo- The transition probability w can be rewritten in 
the following way: 



o PsPO 

27r<5(fc+ - kl^^^^)6{k^ 



fel)x 



The last product of the two (5-functions can be written 
under several convenient forms for the calculation of the 
various integrals. The control parameter k (Eg I45|l rises 
after multiplying the nonlinear operator kernel by the ad- 
vection time Tq [Eq. H17|l ]. Furthermore, assuming unmod- 
ified transverse spectra of the form kj^ with q > 2, when 
we writes the system for 5*+, S" and S", in term of the 
variable y, it can be realized that it is scale invariant. 



Assuming power law solutions for the S's, the coefficients 
of the system are independent of the wave vectors af- 
ter integrating over the angles. Because of the integra- 
tion of the delta- functions over the k's, the system is re- 
duced to a differential system involving three ID-spectra 
depending on a single wavenumber since fc|~ ~ "^ii" 
and fcy ~ ■ Before writing the differential system, we 
approximate the theory specifically for the case y < I, 
where we have already seen that the spatial variation of 
the rms magnetic field is smooth compared to e"^ and we 
fix B at its value Bq at the shock front and introduce the 
amplification factor A = Bq/Boo- We set N — X{y)N^ 
for the three spectra, where we introduce N^ such that 
54fc||)27r(g - 2)P{kJ)-i = knVAooN^iku)/ poVl^. 

The evolution system accounts for the case of substan- 
tial pre-amplification by the non-resonant instability. The 
third-order evolution system reads as 



{X+ + X- 



dX+ 
dy 



{X++X-) 



dx_ 

dy 



{X+ + X-) 



dX" 
dy 



-e~yx+ 


+ X 


[X+X- 4 


-{X+- 


e-yx- - 


Ak. X 


[X+X- 4 


-{X+^ 


Ak, 




2~ 





x[X+ X- + {X+ - X-)X'] (A.: 



For K = 0, the ratio of the first two equations leads to 
X+X- = constant = X+{y = l)X-{y = 1). 

Appendix B: A more complete nonlinear theory 

Because of the efficiency of the backscattering process 
when «; ~ 1, it is reasonable to envisage a secondary gen- 
eration of backward sound waves from backward Alfven 
waves, which also regenerates the forward Alfven spec- 
trum: A~ ^ A+ + S~ . 



N- = - 



d^k^ 



d^k" 



rW 



iV^ 



N"- = 



(27r)3 J (27r)3 

X {N-N+ + N-N"- - N+N"-) 
d^k- f d3fc« 



(27r)3 J (27r)3 

x{N-N+ + N-N 
d^k+ r d^k- 
(2^ 

N-N-' 



w 



(27r)3 

{N-N+ 



w 



N+N-'-) 



N+N"-) . 



(B.l) 



(B.2) 



(B.3) 



We combine the primary and the secondary process, in- 
clude damping of the stable waves (actually the backward 
waves are damped by the cosmic ray streaming at the 
same rate as the forward waves are amplified) . The damp- 



ing rate of the sound waves is y j^k^^Cs- We proceed as 
in the previous appendix to describe the nonlinear evolu- 
tion of the resonant instability, and found that the spectra 
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proportional to S** are still recovered. We then form the 
four order differential system on the amplitudes of the 
spectra that governs the y-profiles: 

dx+ 



{X+ + X- 

{X+ + X-) 
{X+ + X-) 
{X++X-) 



dy 

dX- 
dy 

dX+ 
dy 

dy 



-e-yX+ -Akx 
{X- -X+){Xt + X;) (B.4) 

e~yx^ + Akx 

{X- - X+){Xt + XJ) (B.5) 

{X-X+ + X+X+ - X-X+) (B.6) 
9sX- - ^1 X 

{X-X+ + X'X- - X+X-) (B.7) 



In the case e = 0, is a pure number: Qs = ^| ^"p^ 

which has to be compared with k; typically gs ~ (10~^ — 
10^^)k. It turns out that a relaxation of the soimd waves 
is possible only for > k^{X^ — X~Y, which implies 

X+ ~ X-. 
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